Statistical aspects of nuclear coupling to continuum 
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Various global characteristics of the coupling between the bound and scattering states are ex- 
plicitly studied based on realistic Shell Model Embedded in the Continuum. In particular, such 
characteristics are related to those of the scattering ensemble. It is found that in the region of 
higher density of states the coupling to continuum is largely consistent with the statistical model. 
However, assumption of channel equivalence in the statistical model is, in general, violated. 
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Relating properties of nuclei to the ensembles of ran- 
dom matrices jy is of great interest. A potential agree- 
ment reflects those aspects that are generic and thus 
do not depend on the detailed form of the Hamiltonian 
matrix, while deviations identify certain system-specific, 
non-random properties of the system. On the level of 
bound states the related issues are quite well explored 
and documented in the literature (111- ^ n many cases, 
however, the nuclear states are embedded in the contin- 
uum and the system should be considered as an open 
quantum system. Applicability of the related scatter- 
ing ensemble of non-Hermitian random matrices has 
however never been verified by an explicit calculation due 
to serious difficulties that such an explicit treatment of 
all elements needed involves. These include a proper han- 
dling of multi-exciton internal excitations, an appropriate 
scattering asymptotics of the states in continuum and a 
consistent and realistic coupling among the two. The 
recently developed [[| advanced computational scheme 
termed the Shell Model Embedded in the Continuum 
(SMEC) successfully incorporates such elements and will 
be used below to study conditions under which the sta- 
tistical description of the continuum coupling applies. 

Constructing the full SMEC solution consists of three 
steps. In the first step, one solves the many-body prob- 
lem in the subspace Q of (quasi-)bound states. For that 
one solves the multiconfigurational Shell Model (SM) 
problem : HqqQi = Ei&i , where Hqq = QHQ is the 
SM effective Hamiltonian which is appropriate for the 
SM configuration space used. For the continuum part 
(subspace P) , one solves the coupled channel equations : 
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where index c denotes different channels and Hpp = 
PHP. The superscript (+) means that boundary con- 
ditions for incoming wave in the channel c and outgoing 
scattering waves in all channels are used. The channel 
states are defined by coupling of one nucleon in the scat- 
tering continuum to the many-body SM state in (iV — 1)- 
nucleus. Finally one solves the system of inhomogeneous 
coupled channel equations : 
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with the source term Wi which is primarily given by the 
structure of N - particle SM wave function <!>i and which 
couples the wave function of TV-nucleon localized states 
with (TV — l)-nucleon localized states plus one nucleon 
in the continuum 0. These equations define functions 
lu ■ + ^ , which describe the decay of quasi-bound state $i 
in the continuum. 

The resulting full solution of SMEC equations is then 
expressed as [raj7| : 
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defines the effective Hamiltonian acting in the space of 
quasibound states. Its first term reflects the original di- 
rect mixing while the second term originates from the 
mixing via the coupling to the continuum. Gp is the 
Green function for the single particle (s.p.) motion in the 
P subspace. This external mixing is thus energy depen- 
dent and consists of the principal value integral and the 
residuum : 
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These two terms prescribe the structure of the real W R 
(Hcrmitian) and imaginary W 1 (anti-Hcrmitian) parts of 
W, respectively. The dyadic product form of the second 
term allows to express it as 
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where the M x A matrix V = {Vf} denotes the am- 
plitudes connecting the state (i — 1,...,M) to the 
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reaction channel c (c = 1,...,A) (|].This form of W 1 
constitutes a starting point towards statistical descrip- 
tion of the related effects. In the latter case one assumes 
that the internal dynamics is governed by the Gaussian 
orthogonal ensemble (GOE) of random matrices. Rela- 
tion of this assumption to the classical chaotic scattering 
can also be traced ||. The orthogonal invariance argu- 
ments then imply that the amplitudes Vf can be assumed 
to be Gaussian distributed and the channels independent 
j|. Assuming, as consistent with the statistical ensem- 
ble, the equivalence of the channels one then arrives at 
the following distribution of the off-diagonal matrix ele- 
ments of W 1 for A open channels : 
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with {(wijY 

function. 



r(A/2)^F2( A - 1 )/ 2 
= A. K\ denotes here the modified Bessel 




FIG. 1. Typical projections of distribution of matrix 
W elements coupling to the one channel continuum in 
the SM basis of J 1 = 0+ ,T = states in 24 Mg (his- 
tograms). The projections on imaginary (left) and real axis 
(right) are normalized and plotted versus normalized variable 

a = (W*-(W$))/ox, where o x = {W* 2 ) 1 ' 2 , and X = I,R 
denotes imaginary and real parts respectively. In the upper 
parts all 325 states were taken into account, while in the lower 
parts only 205 states in the middle of the spectrum were in- 
cluded. The full curves represent A = 1 distribution (Eq. 

§)■ 

The physics to be addressed below by making use of 
the above formalism is that of a nucleus decaying by 
the emission of one nucleon. As an example, 24 Mg is 
taken with the inner core of 16 O and the phenomeno- 
logical sd-shell interaction among valence nucleons [ p)| |. 
For the coupling between bound and scattering states a 



combination of Wigner and Bartlett forces is used, with 
the spin-exchange parameter (3 = 0.05 and the overall 
strength coupling = 650MeV • fm 3 || . The radial 
s.p. wave functions in the Q subspace and the scatter- 
ing wave functions in P subspace are generated from the 
average potential of the Woods-Saxon type 

In the above SM space, the 24 Mg nucleus has 325 J 71 " = 
+ ,T = states. Depending on the particle emission 
threshold, these states can couple to a number of open 
channels. Such channels correspond to excited states in 
the neighboring N — 1 nucleus. 
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FIG. 2. Variance of real (o\r) and imaginary (crj) parts of 
matrix elements Wij for one open channel and correlation co- 
efficient p = ((W% Wlj)- (W$) (W^))/{ct r an) between them. 
Different line styles correspond to different daughter nucleus 
spins: 1/2 (full line), 3/2 (dashed line), 5/2 (dot-dash) and 
7/2 (dots). All these quantities are shown as a function of 
energy of the particle in the continuum. 

When testing validity of the statistical model it is in- 
structive to begin with one open channel and to com- 
pare the distribution of the corresponding matrix ele- 
ments with the formula for A = 1. In the example 
shown in Fig. [I], the open channel corresponds to spin 
1/2 and its energy to about the middle of the spectrum. 
Both the imaginary (left) and real (right) parts of W 
are displayed. The upper part of Fig. |l| involves all 325 
J n = 0+,T = states of 24 Mg. Clearly, there are too 
many large and also too many small matrix elements as 
compared to the statistical distribution (solid line) with 
A = 1 . This may originate from the fact that many 
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states in the Q space are localized stronger than allowed 
by the GOE. It is actually natural to suspect that this 
may apply to the states close to both edges of the spec- 
trum. Indeed, by discarding 60 states on both ends of 
the spectrum (205 remain), the picture changes signifi- 
cantly as illustrated in the lower part of Fig. |l|. In this 
case the statistical distribution provides a good repre- 
sentation, interestingly, also for the real part although 
applicability of the formula (0) is not directly justifiable 
as for the imaginary part. Similar behavior is found for 
majority of channels except for a limited number of them 
located at the edges of the spectrum. Hence, the assump- 
tion about the Gaussian distribution of amplitudes Vf is 
well fulfilled in a generic situation. 
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FIG. 3. Fhe same as lower part of Fig. |l| but for 10 chan- 
nels of spins ranging from 1/2 to 7/2 and two energies of the 
particle in the continuum (depicted in the figure). The full 
curves represent V\ e{{ fits with A off indicated while the dashed 
curves correspond to distribution with A = 10. 

As for the equivalence of channels, the conditions are 
expected to be more intricate, especially when differ- 
ent channel quantum numbers are involved, because the 
effective coupling strength depends on those quantum 
numbers. In addition, such a coupling strength depends 
also on energy E of the particle in continuum so the pro- 
portions among the channels may vary with E. This is 
illustrated in Fig. ^| which shows the energy dependence 
of the standard deviations of distributions (as in Fig. [I]) 
of relevant matrix elements for several different channel 
spin values, both for real and imaginary part of W, and 
their correlation coefficient. It should to be noted how- 
ever that within a given spin the differences are much 
smaller. Instead of trying to identify (with the help of 
Fig. ^|) a sequence of A approximately equivalent chan- 
nels and to verify the resulting distribution of matrix 



elements of W against formula (|7]) we find it more infor- 
mative to make a random selection of such channels. An 
example for A = 10 and two different energies (E = 20 
and 40 MeV) of the particle in the continuum is shown in 
Fig. 0. Among these 10 randomly selected channels, two 
correspond to spin 1/2, three to spin 3/2, three to spin 
5/2 and two to spin 7/2. The distributions significantly 
change as compared to those of the lower part of Fig. []]. 
Moreover, Va =w {W^ R ) (Eq. ) (dashed lines) does 
not provide an optimal representation for these explic- 
itly calculated distributions. For E = 20 MeV particle 
energy (the upper part of Fig. |3|), the best fit in terms 
of the formula @) is obtained for A cff = 3.1 for the imag- 
inary part and A cff = 4.4 for the real part of W . At 
E = 40 MeV one obtains A cff = 4.8 and A off = 3.1, cor- 
respondingly. This, first of all, indicates that effectively 
a smaller number of channels is involved what is caused 
by the broadening of the width distribution as a result of 
the non-equivalence of the channels JTT) . Secondly, such 
effective characteristics depend on the energy of particle 
in the continuum, what in turn is natural in view of the 
dependences displayed in Fig. ||. It is also interesting to 
notice that obeys functionally similar distribution as 
W/j although this does not result from Eq. (||) [§]■ 

The fact that generically A off is much smaller than the 
actual number of open physical channels can be antici- 
pated from their obvious non-equivalence in majority of 
combinations as can be concluded from Fig. ||. The 
global distribution, especially in the tails, is dominated 
by stronger channels. Due to the separable form of W, 
which in terms of A explicitly expresses its reduced di- 
mensionality relative to Hqq, an interesting related ef- 
fect in the eigenvalues of Hqq may take place. For a 
sufficiently strong coupling to the continuum one may 
observe a segregation effect among the states, i.e., A of 
them may separate from the remaining M — A states [[l2| . 
This effect is especially transparent when looking at the 
structure of W 1 . For the physical strength °f the 
residual interaction in 24 Mg this effect is negligible, as 
shown in the upper panel of Fig. ^. Only one state in 
this case separates from all others by acquiring a larger 
width. A magnification of the overall strength °f t ne 
coupling to the continuum by a constant factor / allows 
further states to consecutively separate. For / = 7, all 
10 states become unambigously separated as illustrated 
in the middle panel of Fig. ^. Their distance from the re- 
maining, trapped states reflects approximately the order 
of their separation when / is kept increasing. This nicely 
illustrates the degree of non-equivalence of the channels 
and the fact that A off w 5, as consistent with Fig. |^ at 
E — 40 MeV, is an appropriate representation for an ef- 
fective number of relevant open channels. It needs also 
to be noticed that the segregation effect takes place also 
in the direction of the real energy axis, though in this 
sense only three states uniquely separate (again consis- 
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tent with A off = 3.1 of Fig. |3|). This direction of the 
separation originates from the real part of W. Incorpo- 
rating an equivalent multiplication factor into W 1 only, 
results in a picture as shown in the lower panel of Fig. ^. 
No separation in energy can now be observed anymore. 
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FIG. 4. 205 complex eigenvalues for 10 channels and energy 
of the particle in the continuum of 40 MeV are presented as 
small circles with coordinates of Er and Tr. The upper part 
represents those for the original residual interaction between 
Q and P subspaces. The middle one is obtained for 7 times 
stronger interaction , and in the lower part this stronger force 
is applied to W 1 only. 

In summary, the present study indicates that certain 
characteristics of the statistical description of nuclear 
coupling to the continuum, like the distribution of cou- 
pling matrix elements for one channel continuum, do 
indeed apply when the non-generic edge effects are re- 
moved. On the other hand, in realistic SMEC calcu- 
lations we find the generic nonequivalence of channels 
which contradicts the orthogonal invariance arguments 
and results in strong reduction of the number of effec- 
tively involved channels. The quantitative identification 
and understanding of this effect may turn out to be help- 
ful in postulating improved scattering ensembles which 
automatically account for this effect, similarly as vari- 
ous versions of the random matrix ensembles invented 
L3|Jl4|l in the context of bound states. Up to now the 



statistical models ignore the real part of the matrix con- 
necting the bound states to the scattering states. The 
real part of Hqq is likely to be dominated by -Hqq, 
therefore, this, in many cases, may be not a bad ap- 
proximation. Keeping in mind a relatively strong energy 
dependence of W R (see Fig. ||) this however may not be 
true in some cases, especially, because the segregation of 
states in energy (along the real axis) originates from this 
part. Interestingly, W R is found to obey similar statis- 
tical characteristics as W 1 . This does not however yet 
mean that the two parts of W can simply be drawn as 
independent ensembles. In fact, the individual matrix 
elements and W R are often strongly correlated and 
the degree of correlation depends on energy of the par- 
ticle in the continuum. A more detailed account of such 
correlations will be presented elsewhere. 
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